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Abstract. Elements of kinematical and dynamical modeling of ellipti- 
cal galaxies are presented. In projection, NFW models resemble Sersic 
models, but with a very narrow range of shapes (m = 3 ± 1). The 
total density profile of ellipticals cannot be NFW-like because the pre- 
dicted local M/L and aperture velocity dispersion within an effective 
radius (Re) are much lower than observed. Stars must then dominate 
ellipticals out to a few R e . Fitting an NFW model to the total den- 
sity profile of Sersic+NFW (stars+dark matter [DM]) ellipticals results 
in very high concentration parameters, as found by X-ray observers. 
Kinematical modeling of ellipticals assuming an isotropic NFW DM 
model underestimates M/L at the virial radius by a factor of 1.6 to 
2.4, because dissipationless ACDM halos have slightly different density 
profiles and slightly radial velocity anisotropy. In A-body+gas sim- 
ulations of ellipticals as merger remnants of spirals embedded in DM 
halos, the slope of the DM density profile is steeper when the initial 
spiral galaxies are gas-rich. The Hansen & Moore (2006) relation be- 
tween anisotropy and the slope of the density profile breaks down for 
gas and DM, but the stars follow an analogous relation with slightly 
less radial anisotropics for a given density slope. Using kurtosis (hi) 
to infer anisotropy in ellipticals is dangerous, as In is also sensitive to 
small levels of rotation. The stationary Jeans equation provides ac- 
curate masses out to 8-R e - The discrepancy between the modeling of 
Romanowsky et al. (2003), indicating a dearth of DM in ellipticals, and 
the simulations analyzed by Dekel et al. (2005) , which match the spec- 
troscopic observations of ellipticals, is partly due to radial anisotropy 
and to observing oblate ellipticals face-on. However, one of the 15 so- 
lutions to the orbit modeling of Romanowsky et al. is found to have an 
amount and concentration of DM consistent with ACDM predictions. 
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1 Introduction 

The quantity of dark matter lying in the outskirts of luminous elliptical galaxies 
is a hotly debated topic (see Romanowsky, Napolitano, Stoehr, in these proceed- 
ings). There is a wide consensus that, given their flat rotation curves, spiral galax- 
ies must be embedded within dark matter halos, unless one resorts to modifying 
physics (e.g. MOND, see McGaugh in these proceedings). Moreover, dissipation- 
less cosmological TV-body simulatio ns lead to str uctures, whose halos represent 
most (Hayashi et el. 2004) if not all l)Stoehrll200fih spiral galaxies. 

If elliptical galaxies originate from major mergers of sp iral galaxies I Toomrel 



it elliptical galaxies originate irom major mergers oi spiral galaxies ( loomrc 
1977tlMamoilll992tlBaugh et al.lll995 ISprineel et al.ll200l|) . then they too should 
possess dark matter halos. However, inferring the presence of dark matter halos 
in ellipticals is difficult, because the velocity dispersions of the usual kinematical 
tracer, stars, can only be measured out to 2 R e (effective radii, containing half the 
projected light). Moreover, the mass distribution depends on the radial variation 
of the velocity anisotropy, and one cannot solve for both, unless one assumes 
no rotation and makes use of the 4th order moment (kurtosis) of the velocit y 
distribution llvan der Marel fc Franxlll993t iGerhardlllfia lLokas fc Marrionl feooai . 
Using planetary nebu lae (PNe) as tracers of the dark matter at large radii, 



using planetary nebulae (riNej as tracers oi tne dark matter at large radii, 
iRomanowskv et all 1 200.^ found low velocity dispersions for their outermost PNe, 
which after some simple Jeans modeling and more sophisticated orbit modeling 
led them to conclude to a dearth of dark matter in ordinary elliptical galaxies. 
This result is not expected in the standard model of formation of structure in the 
Universe and of ga laxies in particula r. 



This has led us l|Dekel et al.l2005^ to analyze t he final outputs of TV-body simu- 



lation s of spirals m erging into ellipticals 1 Cox et al.l2004l Stoehr, in these proceed- 



ings). iDekel et "all show that the line-of-sight velocity dispersio ns of their simulated 
merger remnants are as low as the PNe dispersions measured bv lRomanowskv et alJ 
thus removing an important thorn in the standard model of galaxy formation in 
a ACDM Universe. One is led to wonder how could the kinematical analysis of 
iRomano wsky et all lead to a lack of dark matter, while the dynamical analysis 
of iDekeietalT nialches the same set of observations with simulations including 
normal amounts of dark matter. 

To explain this discrepancy, it is useful to focus first on some important aspects 
of kinematical modeling. 



2 Kinematical modeling 

2.1 Is the total density profile NFW-like? 

iMerritt et al.l l|2005l) have reported that the projected density profiles of th e halos 
found in dissipationless cosmological TV-body simulations are well fit by ISersid 
( 1968) models, with the shape parameter m»3±l, T his result is not surprising, 
becau se these halos are known to be fairly well fit by the | Navarro. Frenk. fc White! 
(1996, hereafter NFW) model, and lLokas fc Mamonl l(200 1|) had already shown that 
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projected NFW models resemble Sersic models with m = 3 ± 1. 

One should not conclud e that elliptical gal axies, which are known to have Sersic 
surfac e brightness profiles l|Caon et alJflQQsl . are M/L — est NFW models, be- 
cause l|Lokas fc MamorfeOOll) : 1) the range of shape parameters is much narrower 
for the simulated halos (to k, 3 ± 1) than for observed ellipticals (1 < to < 6 
or 10); 2) the normalization of a divergent NFW model with a convergent Sersic 
model leads to ab surd total M/L at i? e ; 3) given the increase of to with luminosity 
llCa.ori et al.ll993l) and the decrease of halo co ncentration with mass l|Navarro et all 
ll997tl,Ting fc Sutoll200(i iBullock et al]l200lh. the trend of increasing best-fit to for 
increasing concentration IjLokas fc Mamonll2f)f)lh implies the unlikely trend that 
elliptical galaxy luminosity decreases with increasing halo mass. 

Now, it is not clear whether the NFW-like structures found in dissipationlcss 
cosmological simulations represent the dark matter component, i.e. are set by the 
dissipationlcss nature of dark matter, or if they represent the total matter, i.e. are 
set by the global gravitational field. However, if one assumes that the NFW-like 
models repre sent the total mass dens ity of elliptical galaxies, one runs into trouble 
in two ways l|Mamon fc Lokasll2005a|k 1) as R decreases below R e , the local M/L 
falls to values well below the M/L of the stellar populations representing ellipti- 
cals; 2) the aperture (R < R e /10) velocity dispersions are ver y much below the 
obser ved values (left plot of Fig. ^| for a given luminosity (the iFaber fc Jackson! 
Il97fil relation). This occurs with th e NFW model, but also with the more recent 
3D-Sersic model that lNavarro et aTl l)2004l hereafter Nav04) found to fit better the 
simulated halos, as well as the steeper mod e ls with inn er slope of —3 / 2 prop osed by 
iFukushige fc Makinol (|l997|) : iMoore et all l|l999|) and ljing fc Sutol (j2000f> . Given 
the limited spatial resolution of cosmological simulations (symbols in Fig. , the 
observed velocity dispersions could be recovered if the total density profile sharply 
steepens (to a slope steeper than —2) just at R e , but this seems a little far fetched. 
The low predicted local M/Ls and aperture velocity dispersions imply that NFW- 
like models cannot represent the total matter. The simplest explanation is that 
these models represent the dark matter component only and that a more concen- 
trated stellar component must dominate NFW-like dark matter within R e . 

The right-hand plot of Figure shows that fitting an NFW model to the total 
density profile of a system composed of Sersic stars and NFW dark m atter yields 
a high concentration, as fit by X-ray observers fe.g-. ISato et aljf2000|) . and the fit 
cannot be very good. 



2.2 Weighing the dark matter out to the virial radius 

The key physical parameter for the dark matter is its mass within the virial radius, 
r v , or cquivalently its M/L. Alas, given that observations are very hard to obtain 
beyond 5 R e (« 0.06 r v ), even with external mass tracers such as PNe and globular 
clusters, o ne requires a strong extra polation to weigh the dark matter component 
out to r„. iMamon fc Lokasl (|2005bh find that the observed velocity dispersion at 
5i? e is a very weak function of M/L at r v . Conversely, M/L at r v is a strong 
(8th power) function of the observed o v at 5 R e . Figure illustrates how badly 
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R/R. 



R/R, 



Fig. 1. Left: Velocity dispersion averaged over a circular aperture of width R for mod- 
els where the total mass density is NFW (dashed), Jing-Suto (2005, their —3/2 slope 
model, dotted) and Navarro et al. (2004, solid) curves, for Lb = Lb,* = 1-88 x IQ 10 L@ 
and Mvir/Ls ~ 39 (lower curves), 390 (middle curves), and 3900 (upper curves) Mq/Lq. 
The shaded region indicates the observational constraints (Faber- Jackson 1976 relation). 
Right: density profiles of an elliptical galaxy made of a Sersic stellar and a c = 9 NFW 
dark matter component. When fitting an NFW to the total density profile (between the 
radii of the two vertical lines), a much higher concentration parameter (c = 35) is found. 



one underestimates M/L at r v by assuming an NFW dark matter model instead 
of the latest 3D-Sersic model of iNavarro et alJ <|2004^ and neglecting radial veloc- 
ity anisotropy . Assuming an isotropic NFW model instead of the slightly radial 
INavarro et alJ model that matches much better the halos in cosmological iV-body 
simulatio ns the M/L derived at r v is 60% too low. But if the orbits are as radial as 
found bv lDekel et al.l in the merger simulations (see Stoehr, in these proceedings), 
then M/L is underestimated by a factor 2.4. 



3 Dynamical modeling 

We have analyzed the end products of -/V-body+gas (SPH) simulations of merging 
equal mass spirals, made of a disk, a bulge, a gas disk, and a dark matter halo 
(see Stoehr et al., in these proceedings). The merger remnants not only show 
stellar surface density profiles that almost perfectly match the observed surface 
brightness profiles of ellipticals (an old result), but more interestingly, they display 
the same line-of-sight stellar velocity dispersion profiles as the st ellar and PN 
velocity dispersion profiles observed in ellipticals l|Dekel et al.ll2005l see Stoehr et 
al., in these proceedings). These simulations highlight how the dark matter and 
stellar kinematics are decoupled. The stellar component is found to dominate 
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Fig. 2. M I Lb at the virial radius divided by the value obtained assuming an NFW dark 
matter model and isotropy, as a function of observed line-of-sight velocity dispersion at 
5 R e (adapted from Mamon & Lokas 2005b). The dark matter model is shown with 
numbers representing the anisotropy radius (where = 1/4) in units of r v . 

the dark matter component within R e whereas one expects from the dark matter 
profiles of cosmological TV-body simulations and the observed s urface brightness 
profile s, that the dark matter only dominates the stars at 2 R e l|Mamon fc Lokasl 
Thus the dark matter has readjusted itself in the inner regions where the 
stars dominate the gravitational potential. Figure |3] shows that for runs with a low 
initial gas content the dark matter slopes are in good agreement with those seen 
in halos of dissipationless ACDM simulations. However, the runs with a high gas 
content show higher slopes, as qualitatively expected since the dark matter should 
respond to the more predominant inner baryons (see Gnedin, these proceedings). 



lHansen &: Moord ((2006) noticed that in dissipationless simulations, (both cos- 
mological and non-cosmological) , the radial profiles of velocity anisotropy and 
density slope are strongly correlated, which should, in principle, allow the lifting 
of the mass/anisotropy degeneracy (see Sec. 1). Figure0]shows that in the merger 
simulations of spiral galaxies made of stars, dark matter and gas, analyzed by 
Dekel et all each component obeys its own correlation (or lack thereof) between 
anisotropy and density slope. The gas relation shows a rotating inner disk, and 
is very dispersed, presumably because of its dissipative nature. The dark matter 
relation is is otropic even in reg ions of steep density profiles, in contrast with the 
prediction of lHansen fc MooreL Interestingly, the stellar relation is close to that 
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Fig. 3. Radial profiles of dark matter density slopes from 10 major merger simulations: 
upper curve & region: gas-rich runs, lower curve & region: gas-poor runs. Dashed 
and solid curves show the predicted slopes for NFW and Navarro et al. (2004) dark 
matter models, respectively, with the concentration parameters given in Mamon & Lokas 
(2005a). 



of lHansen fc MooreL but with over double its dispersion and with (3 roughly 0.15 
lower. Presumably the stellar kinematics are affected by the high dispersion of the 
dissipative gas kinematics from which a few stars are formed. 

The simulations also help to understand what drives the 4th order, gaussian- 
weighted, velocity moment, /14. Figure displays the velocity distribution and /14 
values of the PNe observed in two elliptical galaxies as well as the particle veloc- 
ity distributions in the same simulation run seen in two orthogonal projections. 
Interestingly, the same snapshot viewed along two orthogonal directions produces 
either flat-topped or cuspy velocity distributions. Although the 3D configuration 
is radially anisotropic beyond 1.5 R e , the /14 parameter is either negative (left plot) 
or positive (right plot). The negative /14 is caused by a small amount of rotation 
in the system viewed perpendicular to the angular momentum axis. 

4 Understanding the low outer velocity dispersions 

4. 1 Lack of equilibrium ? 

Although the general Jeans equation in spherical symmetry can be written 



dt dr r 



(W) 



v cot e 5$ 

H V r Vg = -V — , 

r or 
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Fig. 4. Velocity anisotropy vs. logarithmic slope of the density profile in 10 major meger 
simulations. Circles, diamonds and triangles show dark matter, star and gas particles, 
respectively. Filled and open symbols show inner (7? < 2R e ) and outer (R > 2 R e ) 
particles, respectively. Each symbol represents a spherical shell of one of the 10 runs 
analyzed. The shaded regions illustrate /3 = 1 — 1.15 (1 + a/6) ± 0.1 (a is the slope) as 
found by Hansen & Moore (2006) in very different dark matter only simulations. The 
thin horizontal lines represent isotropy. 



when kinematical modelers of observations write the Jeans equation, they always 
omit the terms in v r v$ and the term with a partial time-derivative. With these 
omissions, the derived mass in a shell is 



Afj eans (r) 



G 



dim 



d In ? 



2/3 



(4.1) 



Figure plots the ratio of Mj eans (r) / M (r) . One sees that the Jeans mass is equal 
to the true mass within a few percent out to roughly 8 R e . This means that the 
system is close enough to being stationary that one can safely omit the time- 
derivative in the Jeans equation. Beyond 8R e , the crossing time is long and the 
system has not had time to relax to a stationary state. 



146 



Mass Profiles and Shapes of Cosmological Structures 




(v - v) / a v (v - v) / a v 

Fig. 5. Normalized velocity distributions for two observed galaxies (top) and two snap- 
shots of the same simulated merger remnant (stellar component) viewed perpendicular 
to (lower left) and along (lower right) the angular momentum vector. Short-dashed and 
long-dashed curves represent 0.5 R e < R < 1.5 Re and R > 1.5 i? e , respectively, while the 
solid curve is the gaussian distribution. 

4.2 How much dark matter from orbital modeling? 

While iRomanowskv et all l|2003j) announced that their orbit modeling of NGC 
3379 yielded no solutions with appreciable amounts of dark matter, a closer inspec- 
tion of their 15 orbit solutions (Fig. 0), shows that 2 solutions have the expected 
relatively high M/L b at the virial radius, one of which has the expected dark mat- 
ter concentration, while an additional solution later found by Romanowsky (priv. 
comm.) ha s an even higher M l Lb- Therefore, the kinematic data (stars+PN) 
analyzed bv lRomanowskv et al\ is in fact consistent with the ACDM predictionsl 
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orbit solutions for the elliptical galaxy NGC 3379 of Romanowsky et al. (2003) (circles) 
and an additional solution later found by A. Romanowsky (2005, private communication, 
diamond), adapted from Mamon & Lokas (2005b). For NGC 3379, we assume an absolute 
magnitude Mb = —20.12, derived from a Sersic fit to the B-band surface brightness 
profile of de Vaucouleurs & Capaccioli (1979), and corrected for galactic extinction. The 
vertical shaded region is the recent ACDM constraint on M/L within the virial radius by 
Eke et al. (2005) . The oblique shaded region is the constraint on the ACDM dark matter 
concentration mass relation of Bullock et al. (2001), as rederived by Napolitano et al. 
(2005) for as — 0.9. For both relations, we assume a 40% uncertainty. 
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